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Abstract — An unknown-position sensor can be localized if there 
are three or more anchors making time-of-arrival (TOA) mea- 
surements of a signal from it. However, the location errors can be 
very large due to the fact that some of the measurements are from 
non-line-of-sight (NLOS) paths. In this paper, we propose a semi- 
definite programming (SDP) based node localization algorithm 
in NLOS environment for ultra-wideband (UWB) wireless sensor 
networks. The positions of sensors can be estimated using the 
distance estimates from location-aware anchors as well as other 
sensors. However, in the absence of LOS paths, e.g., in indoor 
networks, the NLOS range estimates can be signiflcantly biased. 
As a result, the NLOS error can remarkably decrease the location 
accuracy. And it is not easy to efficiently distinguish LOS from 
NLOS measurements. In this paper, an algorithm is proposed that 
achieves high location accuracy without the need of identifying 
NLOS and LOS measurement. 

Index Terms — Wireless sensor networks, non-line-of-sight 
(NLOS), time-of-arrival (TOA), semi-definite programming 
(SDP). 



I. Introduction 

Localization algorithms for wireless sensor networks 
(WSNs) have been designed to find sensor location informa- 
tion, which is a major requirement in many applications of 
WSNs. Examples of such applications include animal tracking, 
mapping and location-aided routing. 

Generally speaking, based on the type of information pro- 
vided for localization, protocols can be divided into two 
categories: (i) range-based and (ii) range-free protocols H]. 
Due to the coarse location accuracy, solutions of range-based 
localization are often more preferable and accurate than those 
of range-free schemes. Range estimates from anchors can be 
obtained using received signal strength (RSS), angle-of-arrival 
(AOA) or time-of-arrival (TOA) observations of transmitted 
calibration signals Impulse-based ultra-wideband (UWB) 
is a promising technology where precise ranging can be 
embedded into data communication, due to its robustness 
in dense multipath environments and its ability to provide 
accurate position estimation with low-data-rate communica- 
tion. In this paper, we focus on the investigation of range- 
based localization algorithms for UWB WSNs. One of the 
main challenges for accurate node localization in range-based 



localization algorithms is non-line-of-sight (NLOS) propaga- 
tion which is caused by the obstacles in the direct paths of 
beacon signals. NLOS will result in unreliable localization 
and significantly decrease the location accuracy if its effects 
are not taken into account. This often occurs in an urban or 
indoor environment. Some localization algorithms that cope 
with the existence of NLOS range measurement have been 
proposed [3] [4] [5|, mostly in cellular networks. Roughly 
speaking, there are two categories of approaches to deal with 
the localization problem in the presence of NLOS propagation. 
The first approach identifies LOS and NLOS information and 
discards the NLOS range information for position estimation. 
The second approach uses all NLOS and LOS measurements 
and provides weighting or scaling to reduce the adverse impact 
of NLOS range errors on the accuracy of location estimates, 
they also assume that the NLOS range estimates have been 
identified. 

The number of anchors is typically limited by practical 
considerations. It might be a waste of resources to discard 
NLOS range measurements. To make best use of all range 
measurements, a computationally efficient semi-definite pro- 
gramming (SDP) approach that effectively incorporates both 
LOS and NLOS range information into the estimate of a 
sensor's location is proposed in this paper. We focus on the 
problem of NLOS mitigation, but do not require to accurately 
distinguish between LOS and NLOS range estimates. Given a 
mixture of LOS and NLOS range measurements, our method 
is applicable in both cases without discarding any range 
information. This method is the only SDP based approach to 
reduce the impact of NLOS on node locahzation in WSNs. 
The main advantages of this approach are as follows. 

1) The statistics of the NLOS bias errors are not assumed 
to be known a priori. 

2) No range information is discarded. 

3) NLOS range estimates are not required to be readily 
distinguished from LOS range estimates through channel iden- 
tification. 

In our proposed approach, we assume the following features 
of UWB TOA-based range estimation: the range bias errors 
in NLOS conditions are always positive and significantly 



larger in magnitude than the range-measurement errors in LOS 
conditions. In the next section, we show that the problem 
of node localization, given range information, can be cast 
into a nonlinear programming. We then use SDP relaxation 
techniques and add an additional measurement error to the 
actual range measurements, resulting in a method that suits for 
both LOS and NLOS range estimates to estimate a sensor's 
location. 

The rest of the paper is organized as follows. Section 
II derives the SDP based localization algorithm in NLOS 
environment. In Section III, an extension model was proposed 
to deal with case when the anchor positions are also uncertain. 
In section IV, simulation results are reported. Section V draws 
the conclusion. 

II. NLOS Localization using SDP 

In this section, an SDP based node localization approach is 
proposed. We first introduce the technical preliminaries of this 
algorithm in subsection A and then formulate the localization 
problem as a nonlinear optimization problem in subsection B. 
An extension for our SDP algorithm to the case considering 
the anchor uncertainties is given in Section III. 

A. Background 

The basic setting of this paper is as follows: There are n 
distinct sensors in whose positions are to be determined 
and m anchors whose positions are known a priori. We use 
Xi e R^, i = 1,2, n to denote the sensors and Xj € R^,j = 
n + l,n + 2, ...,n + m to denote the anchors. We use to 
denote the actual distance between anchor and sensor or sensor 
and sensor, i.e., 

r,;j = ||xi -Xjll ,Vi = 1,2, ...,n,j = l,2,...,m + n. (1) 

In practice, we get measurement information for a subset of 
pairs of nodes, which we denote by E. We use Ei to denote 
the measurement information between sensors and anchors and 
E2 to denote the measurement information between sensors 
and sensors. By definition, E — Ei[J £'2. Notice that this 
measurement could be either LOS or NLOS, and since we do 
not distinguish these two measurements, we do not need to 
separate E by the type of measurement. 

In this paper, we assume that the LOS range measm^ement 

is 

where riij ~ N{0,afj) is the measurement error which 
follows a zero-mean Gaussian distribution with standard de- 
viation (Ti J. 

Similarly, the NLOS range measurement is assumed to be 

A J rij + + Si J, (3) 

where 6ij is the error of NLOS measurement. 

The idea of our approach is to get an upper bound as well as 
a lower bound for the true distance of each pair of nodes (could 
be either anchor and sensor or sensor and sensor) based on the 
measurement we observe, without distinguishing whether it 
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Fig. 1. An instance with two anchors 

comes from LOS or NLOS measurement. These bounds will 
form a feasible region for possible locations of each sensor 
and we then choose one "center" point from this region as 
our estimation. 

First we show how we obtain the upper bound for the 
distance of certain pair of nodes. Since in the NLOS case, 
the measured distance is larger than the actual distance, the 
measurement itself is an upper bound. For the LOS case, we 
have the upper bound as: 

where nfj is an upper bound on the measurement error, 
which could be calculated in advance based on experimental 
measurements. Therefore, we have a uniform upper bound for 
the distance between each pair of nodes as follows: 

r,,, - ||x, - Xj II < d,- , + riY^j V(^, j) G E. (4) 

Next we derive the lower bounds. Here we use the same 
idea as it is in f6|. We first consider the distances between 
sensors and anchors. We consider those sensors that have more 
than two anchors in range. For each of those sensors, we 
draw a circle for each anchor in range centered at the anchor 
position and the radii is the upper bound computed by the 
above method. Obviously, these circles will have a common 
intersection part which contains the true location of the sensor 
Now we look at each pair of the circles. As shown in Fig. 1, 
a lower bound of the distance between the sensor and anchor 
1 is \\di — AB\\, where AB is the intersection of the line 
connecting the centers of the two circles and the common area 
of the two circles. Then we take the maximum of this bound 
over all the anchors in range, and get a final lower bound as 
follows: 

> k,j = niax \\dj - ABkW 

k,{i,k)£Ei 

for any pair of sensor and anchor in range. In the above 
formulation dj is the radii for circle j and ABk is the line 



segment based on the intersection of circles j and k. Therefore, 
for each pair of sensor and anchor in range, we get a lower 
bound for their distances. 

However, for the distances between sensors and sensors, 
we cannot apply the same technique because the position of 
the sensors are not known. Thus for those sensors, the lower 
bound is set to zero. 

Therefore, we get the following upper and lower bounds: 
For each in-range sensors and anchors, we have 

and for each in-range sensors and sensors, we have 

< rij < di,j + 

For the later convenience, we uniformly write the upper and 
lower bound for the distance between node i and j by utj and 
li^j, respectively. 

Remark 1: In some circumstances, we do not have the 
communication between sensors and sensors. In that case, we 
simply remove the constraint between sensors, only keeping 
those between sensors and anchors. 

Remark 2: This approach can also be applied to the cases 
when we have prior information on which measurement is 
from LOS path and which is from NLOS path. If we know 
a priori that a certain measurement is from NLOS path, 
then we can compute the upper bound by simply using 
the measurement, or if we know the error is in a certain 
distribution, then we can again adjust the upper and lower 
bound accordingly. The same thing applies when we know a 
certain measurement is from LOS path. 

B. Localization Algorithm using Semi-definite Programming 
for both NLOS and LOS Environments 

In this section, we present a convex optimization algorithm 
for node localization based on the bounds we obtained in the 
previous subsection. 

As shown in Fig. 2, for the single constraint case (s < 
||x — a|| ^ R), it is easy to see that one heuristic position 
estimate lies on the circle with center a and radius ^^-^ . (By 
way of example, the square indicates the possible position for 
an efficient position estimate in Fig. 2.) 

This can be determined by minimizing the following ex- 
pression: 

(||x-a||-s)2 + (||x-a||-i?)2. (5) 

On expanding, Q becomes 

(||x-a||-.f + (||x-a||-i?)2 
= 2||x - a|p - 2(s + i?)||x - all + + i?2 

where and are constants defined in the previous sub- 
section. 

Therefore, the optimization problem for locating the sensors 
can be formulated as: 

mina; T,^<j■.{^.j)eE[\\^^ - ^jf - '^ikj + Ut.j) ||x, - x^H]. 

(7) 




Fig. 2. The single constraint 

Obviously, (|7]i is nonconvex, which cannot be solved easily. 
However, we can relax the problem to two convex optimization 
problems by using the SDP relaxation techniques as proposed 
in |7| and |8|, which are refeiTed as FullSDP and ESDP, 
respectively. 

Define 

X=[xi,X2,...,X„+™]Gi?2x("+") 

Y = X^X. 
We also define: 

2 

g,j = ||x, -Xjll . 
Then we can write (O as follows: 

s.t. 7i J = Yii + Yjj — 2Yij ^jQ^ 

Y = X'^X h 0. 

By performing the SDP relaxation, we relax (fTOb to a convex 
program as follows: 

S.t. T^ij ^ii ^■^J 

i..i>sl., (11) 

where I2 is the 2x2 identity matrix. 

It is worth noting that the anchor part of Z is known, 
therefore it is also in the constraint. Hence, we formulate the 
localization program in (fTTb . 

When the problem is large, the SDP formulation might be 
slow Based on the work of fS], we can further relax it 



into an ESDP formulation: 



(12) 



where ^(i_2,ij) denotes the principal submatrix of j consisted 
of row and column l,2,i,j. 

Both ( fTTl ) and ( fT2l can be solved by standard SDP solvers 
such as SeDuMi or SDPT3 in a centralized way lllOl . We 
choose YALMIP 1 11 1 as the programming interface. 

Remark 3: In practice, one might want to add different 
weights to different terms in objective according to the confi- 
dence he has on each measurement. For instance, we can give 
a lower weight to NLOS part if we have prior statistics on 
NLOS measurements. 

In the next section, we are going to discuss one extension of 
above model, i.e., the situation in which there are uncertainties 
in the anchor positions. We show that a similar SDP model 
can be formulated to solve this problem. 

III. Localization in NLOS Environment 

CONSIDERING ANCHOR ERRORS 

In the model where uncertainties exist for anchor positions, 
we assume that the true anchor positions are within a certain 
ball around the estimated ones, namely, for each anchor j 



where Xj is the estimated value while Xj is the true one. In our 
cases, Uj is given, and it usually comes from the confidence 
in the measurement uses. 

In this case, in addition to the objective in (|7]l, we also add 
the anchor error in the objective. And the anchor positions 
become variables as well. 

We formulate the problem with anchor position uncertainty 
as follows: 



i<j:{i.,j)£E iW-^^ 

X7" X.> 



2ik 



(13) 

where the second term is about the anchor uncertainty. Note 
that one can also add some weight to each term, denoting the 
different confidence level one has for each measurement. 

By using the same technique, we can relax (T3[ into a 
convex program. Again, define 



X = [xi,X2, ■■■,Xn+: 

h X 

x^ x^x 



,] e i?2x(n+m) 



z = 



and 



9i.j 



llx, 



(14) 



(15) 



Then we can write ( fTsT i as follows: 



EJi + Jn r ry 



S.t. 



h 



(16) 



Z = 



X 



X' x'x 



>- 0. 



By performing the SDP relaxation, we relax (fT6b to a convex 
program as follows: 

s.t. 7i J- — Zi^2,i+2 + ■^j+2 j"+2 — 2Zi+2 J-+2 (-jy-j 

Z{i,2) = h 

By using the same method as in (fT2] i. we also get the ESDP 
relaxation to this case: 

min E»<j;(jj)G£;[7»j - '^ikj + Uij)gij] 

+ jyj=n+l[Zj+2,j+2 - 22'J'2,j+2^j] 
S-t. 7ij- = ^i+2,i+2 + ■^j+2j"+2 ~ 2Zi_|_2_j+2 ^g-j 

Z(i,2) = h 

Zil,2,^,J)hOV{i,J)eE. 

IV. Numerical Results 

In this section, simulation results are presented and an- 
alyzed. The performance evaluation focuses on the posi- 
tion estimation accuracy of the proposed algorithm. We 
consider a 2-dimensional region with a size of 40 m x 
40 m. There are totally 18 anchors locating in the area. 
Eight of them are located at the boundary (20, 20)m, (- 
20,20)m, (20,-20)m, (-20,-20)m, (0,0)m, (-20,0)m and (0,- 
20)m, while the remaining ten anchors are randomly deployed 
in the area. In this simulation they are localized at (4.3416,- 
19.3696)m, (-19.3458,-12.3970)m, (3.4767,-17.6967)m, (- 
5.2972,5.2580)m, (8.7053,7.7067)m, (-16.6368,-1.8257)m, (- 
2.3268,-5.8699)m, (-13.8557,7.0257)m, (7.9685,9. 1003)m and 
(-0.8646,2. 1936)m. Then, we deploy 80 sensors in the field, 
but the number of sensor is too large to be listed, so we omit 
the listing of their coordinates. Nevertheless, the configuration 
is shown in Fig. [3] For the sake of simplicity, we assume that 
all the sensors and anchors can see each others, i.e., the fully 
connected situation is considered. 

We follow the noise model of (|3]l with Ui a normally 
distributed variable with noise power -40dB and Si being a 
uniformly distributed random variable drawn from [0,0.5]. 
That means all measurements contains an NLOS error Our 
proposed SDR formulation of ( fTTT i is applied to find the 
estimated sensor positions and the result is shown in Fig. |4] 
The average mean square position error is 6.66 13m^. From 
the figure we see that the proposed method can provide good 
estimation by mitigating the effects of NLOS measurements. 
The solution can act as an initial guess for other numerical 
search to obtain better estimate. 
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Fig. 4. The configuration of the network 



V. Conclusions 

A semi-definite programming based node localization al- 
gorithm in NLOS environments for UWB wireless sensor 
networks has been proposed in this paper. The problem of node 
localization in the presence of anchor position uncertainty has 
been approximated by a convex optimization problem using 
the SDP relaxation technique. Given a mixture of LOS and 
NLOS range measurements, our method is applicable in both 
cases without discarding any range information. Simulation 
results demonstrate the effectiveness of our method. 
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